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Abstract 

Within the Standard Model, a review of the current determination of the sides and 
angles of the CKM unitarity triangle is presented, using experimental constraints 
from the measurements of \ek \ , \ V u b jV c b \ , Afid and from the limit on Am s , available 
in September 2000. Results from the experimental search for B® — B® oscillations are 
introduced in the present analysis using the likelihood. Special attention is devoted 
to the determination of the theoretical uncertainties. The purpose of the analysis 
is to infer regions where the parameters of interest lie with given probabilities. The 
BaBar u 95% C.L. scanning" method is also commented. 



To be submitted to JHEP 



1 Introduction 



In the Standard Model, weak interactions of quarks are governed by the four parameters 
of the CKM matrix [[3J which, in the Wolfenstein parametrisation ||, are labelled as 
A, A, p and fj []. Measurements of semileptonic decays of strange and beauty particles are 
the main sources of information on A and A, respectively. The values of \sk\, \V u b/V c b\, 
Amd and Am s provide a set of four constraints for p and fj. These constraints depend, in 
addition, on other quantities obtained from measurements and/or theoretical calculations. 
The regions of p and fj preferred by the four constraints are expected to overlap, as long as 
the Standard Model gives an overall description of the various experimental observations. 

Since several years there has been intense activity to constrain the allowed region in 
the (p, fj) plane from the best knowledge of the experimental and theoretical inputs [|]]- 
|20fl . Though the analysis methods differ in some details, they have a common ground in 
what we shall call standard approach through this paper. First, the goal of the various 
authors has been, explicitly or implicitly, to infer regions in which the values of p and fj are 
contained with a certain level of probability (or confidence) . Second, uncertainties due to 
statistical errors and systematic effects in experiments, as well as theoretical uncertainties, 
are combined together to deduce a global uncertainty about p and fj. As far as this second 
point is concerned, the various authors have used different "prescriptions" which can be 
seen, indeed, as approximations of the consistent Bayesian method which is described, 
and adopted, in this paper. The 68% probability regions favoured by the data and the 
theoretical understanding of the relevant processes select quite narrow regions for p and 
fj, largely independent of the details of the specific methods and of the different treatment 
of (experimental) systematic and theoretical uncertainties. 

A different approach, named u 95% C.L. scanning in this paperQ, has been adopted in 
the BaBar Physics Book [p3j , and recently used in |24[ . In this approach, it is stated that 
it is not possible to define probability distributions for theoretical parameters coming from 
calculations affected by systematic uncertainties or based on educated guesses (in practice 
all theoretical parameters and some experimental systematics belong to this class). On the 
basis of these considerations, the 95% C.L. scanning approach rejects the two basic points 
of the standard method and a different procedure is proposed. For the theoretical inputs, it 
is assumed that one can only define intervals inside which the true values of the parameters 
are contained. At fixed values of the theoretical inputs (within the allowed intervals) a 
maximum likelihood fit, which includes the other sources of uncertainties, is made, and 
95% C.L. contours are determined. Finally, the envelope of such contours is ^proposed to be 
a (conservative) method to obtain some 95% C.L. regions for all CKM parameters" J24 . 



Using the same arguments of the BaBar Physics Book p3| |, this procedure has been 



recently recommended in |25|, |26] , as opposed to the standard approach which is claimed 
of being too optimistic. 

In view of the importance of constraining the parameters of the CKM matrix, or of the 
possibility of detecting signals of new physics in low-energy weak decays, in this paper 
we reconsider the whole matter, and in particular we focus on the most critical issues 
of the CKM-triangle analysis, namely the uncertainty on theoretical parameters and the 
inferential framework to handle consistently all uncertainties. This also allows to answer 



1 P = P(l-ir) and fj = r)(l-^) 



2 Note that this 95% C.L. scanning is different from the "scanning method" used to predict e.g. e'/e [^jj 
(see aiso [E2J for comments). 
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to several important controversial questions raised in the past, namely: 



• whether it is possible, or necessary, to assign a probability distribution function 
(p.d.f.) to theoretical parameters; 

• whether it is possible to define a p.d.f. for quantities extracted from physical mea- 
surements and from theoretical parameters affected by systematic uncertainties; 

• whether average values, errors and p.d.f. of the quantities considered in the present 
analysis depend in a crucial way on the assumptions made on the theoretical pa- 
rameters and the corresponding errors. 

The main conclusions of this study are the following: 

1. The standard method is a theoretically sound approach, which finds its justification 
within the inferential framework discussed in this paper. This method allows a 
consistent treatment of the systematic and theoretical uncertainties and makes it 
possible to define regions where the values of p and fj (as well as of other quantities 
of phenomenological interest) are contained with any given level of confidence. In 
this respect the criticisms of []25|, are not justified. 



2. The BaBar 95% C.L. scanning, instead, is based on an ad hoc prescription intended 
to define only a "95% C.L." region. The meaning of this statement is unclear: the 
so called "95% C.L. region" does not correspond to the usual statistical definition 
of 95% confidence that the values of parameters lie in that region, neither in a 
frequentist sense, nor in a Bayesian one. 

3. For the sake of comparison, with the standard method we have used the same values 
of the input parameters and tried to mimic the same uncertainties of the 95% C.L. 
scanning. We find that the "95% C.L." regions selected in the (p, fj) plane with the 
two methods are very similar, and it is thus not founded to qualify as too optimistic 
the standard approach. 

4. In the 95% C.L. scanning approach the information contained in the p.d.f. of the 
relevant quantities (p, fj, sin (2/3), etc.) is missing. Thus we only know that a certain 
quantity is somehow expected in a given interval, but we do not know which is the 
most probable value, what is the shape of its p.d.f. etc. We show, instead, that this 
important information can be extracted from the data, using the standard method, 
in spite of the uncertainties in the theoretical parameters. 

Regarding the analysis, since this kind of studies have been extensively illustrated in 
previous publications, see for example fl4| and |L9[ , here we only discuss in detail — B® 



mixing, for which we adopted a different procedure, based directly on the likelihood 
(Section |]). 

The main results for the physical quantities (p, fj, sin(2/3), etc.) can be found in 
Sections 0. 

The remainder of the paper is organised as follows. In Section |2] we summarise the 
theoretical constraints between p and f] and the available experimental and theoretical 
inputs in the Standard Model. In Section |3] we describe the inferential framework used 
in this study and relate it to different standard approach analyses. Our comments on 
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the 95% CL. scanning method are also presented in this Section. The choice of values 
and uncertainties for the most critical theoretical parameters is discussed in Section |J In 
Section [|the p.d.f. determination for \ V c b\, \ V u b\ and the parameter A are explained. A new 
method to include the information coming from searches of B® — B° s mixing is illustrated 
in Section ||. The results of the analysis are presented and discussed in Section |7]. The 
stability of the results has been verified by varying the different input parameters in 
Section ||. A comparison of our results with those obtained with 95% CL. scanning is 
made in Section 0. Finally, conclusions are drawn in Section [TO. 



2 Standard Model formulae relating p and 77 to ex- 
perimental and theoretical inputs 

Four measurements restrict, at present, the possible range of variations of the p and fj 
parameters: 

• The relative rate of charmed and charmless 6-hadron semileptonic decays which 
allows to measure the ratio 



V. 



V F 



cb 



The B% — B® time oscillation period which can be related to the mass difference 
between the light and heavy mass eigenstates of the B® — B® system 

Am d = ^m 2 w r) c S(x t ) A 2 X 6 [(1 - p) 2 + fj 2 ] m Bd f 2 B Bd , (2) 



where S(x t ) is the Inami-Lim function p7]| and x t = m 2 /M^ v . m t is the MS top 
mass, mf' IS (mf IS ), and r] c is the perturbative QCD short-distance NLO correction. 
The remaining factor, f Bd BB d , encodes the information of non-perturbative QCD. 



Apart for p and fj, the most uncertain parameter in this expression is fB d \jBB d - The 
value of i] c = 0.55 ±0.01 has been obtained in p8| and we used m t = (167 ±5) GeV, 
as deduced from measurements of the mass by CDF and DO Collaborations 



The limit on the lower value for the time oscillation period of the B° s — B° s system 
is transformed into a limit on Am s and compared with Am^ 



■ufBa^Ba ( A \ f , , 2 _ 2l , , 

[0--P) +V ]• (3) 



Am d rn Bd f 2 B d B Bd ( A 
Am s " m B J 2 Bs B Bs [l - $ 



The ratio £ = fs s \J B B J fs d \J k>B d is expected to be better determined from theory 
than the individual quantities entering into its expression. In our analysis, we 
accounted for the correlation due to the appearance of Am d in both Equations (|2]) 
and (D. 
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CP violation in the kaon system which is expressed by \ek\ 



\ £ k\ 



where 



C e A 2 \ 6 fj 



-T]iS(x c ) +r] 2 S(x t ) (A 2 A 4 (1 - p)J +r] 3 S(x c ,x t ) 



B 



K 



a 



G]pf 2 Krn K m 2 w 



(4) 



(5) 



S(xi) and S(xi, OC j^j £1X6 the appropriate Inami-Lim functions [|27| of x q = m; l/mly, 



including the next-to-leading order QCD corrections |28 
parameter is B^- 



301. The most uncertain 



Constraints are obtained by comparing present measurements with theoretical expecta- 
tions using the expressions given above and taking into account the different sources of 
uncertainties. In addition to p and -fj, these expressions depend on other quantities which 
have been listed in Table [| Additional measurements or theoretical determinations have 
been used to provide information on the values of these parameters. 



3 Inferential framework 

In this Section we recall the basic ingredients of the standard method, interpreted in the 
framework of the Bayesian approach. This allows us to discuss the role of the systematic 
and theoretical uncertainties in deriving probability intervals for the relevant parameters. 

3.1 Standard approach and Bayesian inference 

Each of Equations (0)-© relates a constraint Cj (where Cj stands for IV^/V^I, Am^, 
Am,j/Am s and \ek\, for j = 1, ... ,4) to the CKM-triangle parameters p and fj, via the 
set of ancillary parameters x, where x = {x\, x 2 , ■ ■ ■ , x^} stand for all experimentally 
determined or theoretically calculated quantities from which the various Cj depend 

C J = Cj(p,fj]x.). (6) 

In an ideal case of exact knowledge of Cj and x, each of the constraints provides a curve in 
the (p, fj) plane. In such a case, there would be no reason to favour any of the points on the 
curve, unless we have some further information or physical prejudice, which might exclude 
points outside a determined physical region, or, in general, assign different weights to 
different points. In a realistic case, we suffer from several uncertainties on the quantities Cj 
and x. Uncertainty does not imply, however, that we are absolutely ignorant about a given 
quantity. First of all, there are values which, to the best of our knowledge, we consider 
ruled out (for example a value of m t of 100 GeV or 500 GeV). Second, we assign different 
probabilities to the values within the "almost certain range", 147 GeV < m t < 187 GeV 
say 0. In the m t case, for example, we think that it is much more probable that the value 
of m t lies between 157 and 177 GeV rather than in the rest of the interval, in spite of the 
fact that the two sub-intervals have the same widths. 

This means that, instead of a single curve (Q) in the (p, fj) plane, we have a family of 
curves which depends on the distribution of the set {cj,x}. As a result, the points in the 

3 In this example m t is the MS top mass of Equation (0), mt = (167 ± 5) GeV. 
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ArriK 


0.5301 x IO" 2 ps" 1 


fixed 
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G F 


1.16639 x 10~ 5 GeV~ 2 


fixed 
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mw 


80.42 GeV 


fixed 
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5.2792 GeV 


fixed 
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5.3693 GeV 


fixed 
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m K 


0.493677 GeV 


fixed 




m 



Table 1: Values of the quantities entering into the expressions of \ek\, \V u b/Vcb\, Amj and 
Am s . In the third and fourth columns the Gaussian and the flat part of the uncertainty are 
given, respectively. 

(p, ff) plane get different weights (even if they were taken to be equally probable a priori) 
and our confidence on the values of p and f] clusters in a region of the plane. 

The above arguments, which we consider very natural and close to physicist intuition, 
can be formalised by using the so called Bayesian approach (see |54| for an introduction). 
In this approach, the uncertainty is described in terms of a probability density function 
/(•), which quantifies our confidence on the values of a given quantity. The inference of p 
and f\ becomes then a straightforward application of probability theory, getting rid of all 
"ad hoc prescriptions". 

The simplest way to implement the probabilistic reasoning discussed above is to define 
an idealised "p.d.f." for each constraint 

/(Ml c i' x ) K S ( c j ~ c i(M x )) , (7) 

where 5 is the Dirac delta distribution. The quotes recall us that this p.d.f. is a distri- 
bution in a mathematical sense, which is to be taken as the limit of a very narrow p.d.f. 
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with values different from zero only along a curve. The p.d.f. which takes into account 
the full uncertainty about Cj and x is obtained from (|7|) by making use of the standard 
probability rules 



f(p,v) = J f(P,v\cj,*) ■ /(cj,x)dcjdx 

oc J 6(cj - Cj{p,f},x)) ■ f(cj) ■ / (x) dcj dx 



'OC 



oc 



S(Cj -Cj(p,T],X)) 



1 



2n a(d 



exp 



[Cj Cj y 



(9) 

•/(x)dcjdx (10) 



1 



27r a(cj 



exp 



2 a 2 '(a) 



2^( c .) 

/(xi) • /(rr 2 ) • • • f(x N ) dx , (11) 



where Cj is the experimental best estimate of Cj, with uncertainty o~{cj). A Gaussian 
distribution has been assumed just for simplicity and without lack of generality. The 
joint p.d.f. /(cj,x) has been splitted as a product of the individual p.d.f., assuming the 
independence of the different quantities, which is a very good approximation for the case 
under study. 

As alternative procedure, one may introduce a global inference relating p, fj, Cj and x. 
This is followed by a second step where marginalization is performed over those quantities 
which we are not interested to. In this case, by making use of Bayes' theorem, we obtain 



f(P,V,Cj,X.\ Cj) OC f{Cj\Cj,p,fi,XL) ■ f(cj,p,fj,x) 

°c f(Cj I Cj) ■ f(Cj I p, fj, x) • /(x, p, fj) 

OC f(Cj I Cj) ■ 5{Cj - Cj{p, fj, x)) • /(x) • f (p, fj) 



(12) 
(13) 
(14) 



where fo(p,fj) denotes the prior distribution as often used in the literature. The var- 
ious steps follow from probability rules, by assuming the independence of the different 
quantities and by noting that dj depends on (p,fj,y:) only via Cj. This is true since Cj is 
univocally determined, within the Standard Model, from the values of p, fj and x (hence 
the limit to a delta function of its p.d.f.). We then recover Equation i) by assuming a 
Gaussian error function for Cj around Cj\ ii) by considering the various Xi as independent; 
iii) by taking a flat a priori distribution for p and fj and iv) by integrating Equation ([14]) 
over Cj and x. Note that equiprobability of all points in the (p, fj) plane was also implicit 
in Equation as discussed above. 

Although the first derivation of Equation (|S|) is probably the most intuitive one, here- 
after we use the second one, which is the usual way of performing Bayesian inference. 
The second procedure also shows explicitly the connection with the methods which we 
denoted as standard in the introduction. 

The extension of the formalism to several constraints is straightforward. We can 
rewrite Equation ([12]) as 



f(p, V, x I ci, 



cm 



^ n fj( c j\p^^) x n fifa) x f°(p,v) 



(15) 



=1,M 



=1,N 



In the derivation of fllq) , we have used the independence of the different quantities. 
Moreover, the conditioning from fj(-) on the Cj have been removed, since the Cj act as 
intermediate variables which are finally integrated away. The derivation (^)-(jT^) can also 
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be easily extended to the case of several constraints and leads, again, to the same result 
as that found by using Bayes' theorem. We have only to account properly the weight on 
Cj, induced by the other constraint (s) previously considered. In this case Equation (|7|) 
becomes f(p, fj \ c x , ■ ■ ■ , c M , x) oc 6(c M - c M (p, 77, x)) • f(p, r] | c x , . . . , c M _ 1; x). 

By integrating Equation (|15D over x we can rewrite the inferential scheme in the 
following convenient way 

f(p, fj | c, f) a £(c | p, fj, f) x f (p, fj) , (16) 
where c stands for the set of measured constraints, and 

£{c\p,fj,f)= ]J fjicj |p,?7,x) n fi(xi)dxi (17) 

J j=l,M i=l,N 

is the effective overall likelihood which takes into account all possible values of Xj, properly 
weighted. We have written explicitly that the overall likelihood depends on the best 
knowledge of all Xi, described by /(x). 

Whereas a priori all values for p and fj are considered equally likely, a posteriori the 
probability clusters around the point which maximises the likelihood. This is the reason 
why, in principle, different procedures for determining p and fj, based on the maximum 
likelihood, are equivalent to the method described here and should get similar results. 
We say "in principle" because other methods are typically implemented using the \ 2 
minimisation. This implies the assumption of a multi-Gaussian solution of the integral 
(|17|), with overall standard deviations which are simply a quadratic combinations of the 
"uncertainties" related to each x^. On the other hand, a quadratic combination relies on 
the approximative linear dependence of Cj from the possible variations of x^. 

In conclusion, the final (unnormalised) p.d.f. obtained starting from a flat distribution 
of p and fj is 

f{p,fj) oc / J] fji&j \p,fj,x) J] fiix^dxi. (18) 

J j=l,M i=l,N 

The integration can be done by Monte Carlo methods, the normalisation is trivial, and 
all moments can be calculated in a (conceptually) easy way. Obviously there are several 
ways to implement the Monte Carlo integration, using different techniques to generate 
events. A comparison of the results obtained with the approach of [O, where some effort 



has been done to improve the generation efficiency, and the results of |14[ is presented in 
Section |7j. 

It is important to note that the inferential method does not make any distinction 
on whether the individual likelihood associated to some constraint is different from zero 
only in a narrow region (and we usually refer to this case as "measurement"), or it goes 
to zero only on one of the two sides (e.g. when Cj — > oo or 0). In the latter case, the 
data only provide an upper/lower bound to the value of the constraint. This is precisely 
what happens, at present, with Am s . Therefore, the experimental information about this 
constraint enters naturally in the analysis (more details can be found in Section |6|). 

3.2 Treatment of systematic and theoretical uncertainties 

At this point, it is in order a discussion on some important ingredients of the analysis 
which raised some controversy in the past. They are related to the quantitative handling 
of the uncertainties due to systematic effects and to theoretical inputs. 
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In Equation (|9[) we have written explicitly that Cj is Gaussian distributed around Cj . As 
a consequence, we tend to say that also Cj is Gaussian distributed around Cj, the inversion 
being well understood in the case of random errors. The question is how to include the 
case where also systematic uncertainties are present. One of the nice features of the 
Bayesian approach is that the uncertainty has positively the same meaning, and there 
is no conceptual distinction between the uncertainty due to random fluctuations, which 
might have occurred in the measuring process, the uncertainty about the parameters of the 
theory, and the uncertainty about influence quantities (i.e. "systematics" ) of not-exactly- 
known value, (see |34|| and f35j). Under the assumption that the individual likelihoods 



fj(-) of Equations ([T5|)-(|l8|) do only depend on random effects, the uncertainty due to 
systematics can be included using, again, concepts and formulae of conditional probability. 
In fact, calling h the set of influence quantities on which the measured constraints may 
depend, with joint p.d.f. g(h), the likelihood ( jT7|) becomes 

£(c|p,77,f,g) = / Yl fj(cj \p,f},x,h) ■ /(x) -#(h)dxdh, (19) 

j=l,M 

where we have written the p.d.f. of x in its general form, allowing also correlations among 
the elements. As can be seen from Equation ([19]), there is neither a conceptual nor a 
formal distinction between the handling of x and h. Therefore, we can simply extend 
the notation to include in x the influence parameters responsible of the systematic uncer- 
tainty, and use Equation (|T7|) in an extended way. This is what it has been actually done 
in the past to infer f(p,f]) with the Monte Carlo integration method resulting from ([18]). 
Moreover, see also and references therein, we arrive to the following conclusion. Ir- 
respectively of the assumptions made on the p.d.f. of x, the overall likelihoods f(£j) are 
approximately Gaussian because of a mechanism similar to the central limit theorem (i.e. 
just a matter of combinatorics). This makes the results largely stable against variations 
within reasonable choices of models and parameters used to describe the uncertainties due 
to theory and systematics. This also explains why methods based on x 2 minimisation (for 



example refs. [[RJ, []T7| and |I8]| ) can be considered as approximations of the one used in 
refs. ||13||,|[14j|, [|15|1 and ]19[. We stress again that a common ground of all the methods 
that we classify as standard is to produce regions where p and f] are contained with any 
given level of confidence. On the contrary, the BaBar 95%C.L. scanning is based on an 
ad hoc prescription which obscures the meaning of the results. In that approach, a so 
called "95% C.L." is produced, which does not correspond to the usual 95% confidence 
that the parameters lie in that regions. 

As far as the choice of the mathematical expression for fi(xi) is concerned, practical 
examples of simple models can be found in |]35 ]. Due to the insensitivity of the result on 



the precise model, as discussed previously and as shown later in this paper, we simplify 
the problem, by reducing the choice only to two possibilities. We choose a Gaussian model 
when the uncertainty is dominated by statistical effects, or there are many comparable 
contributions to the systematics error, so that the central limit theorem applies. We 
choose a uniform p.d.f. if the parameter is believed to be (almost) certainly in a given 
interval, and the points inside this interval are considered equally probable. 

A final comment, before ending this Section, concerns the compatibility among infer- 
ences provided by individual constraints. In the simplified approach based on x 2 min- 
imisation, a conventional evaluation of compatibility stems automatically from the value 
of the x 2 & t its minimum. This information is lost in the likelihood approach, but we 
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accept this loss without any regret, in view of what we gain. It is well known, indeed, 
that crude arguments about compatibility or incompatibility, based only on the minimum 
value of the \ 2 an d the number of degrees of freedom lead to misleading conclusions (as 
premature claims of new physics have demonstrated in the past decades). Therefore, we 
consider more reasonable to judge the compatibility of the constraints by comparing par- 
tial inferences obtained when removing each constraint at the time. Examples are given 



in Figures |9] and [11]. In our analysis, the overlap of the various constraints is excellent, 
and therefore we have no reason to suspect deviations from the Standard Model, given 
the available experimental information. 



4 Theoretical inputs: B K , fB\Bs and the £ param- 
eter 

In this Section we discuss the theoretical inputs which have been used in our phenomeno- 
logical analysis. In particular, we explain how central values and uncertainty models for 
the different quantities have been chosen. This gives us the opportunity of clarifying some 
issues on which there is, we think, some confusion in the literature. This discussion may 
be instructive especially for non-lattice experts (often experimentalists) who are engaged 
in this kind of analyses and have to find some orientation in using results from a plethora 
of lattice studies. 



For example, in ref. |20[| , the renormalization-scheme dependence of the meson de- 
cay constants is mentioned. Unfortunately decay constants (which are related to matrix 
elements of the weak axial current), as all measurable physical quantities, are scheme 
independent by definition. Indeed the physical, scheme-independent axial current is ob- 
tained from the lattice one by a finite renormalization constant, An = ZaA 1 °; u , which 



can be (has been) determined non-pertubatively with a negligible uncertainty |3B|, |37] . 



Scheme dependence only enters some theoretical predictions (but never in those for the 
decay constants) because the perturbative calculations of the Wilson coefficients in the 
effective Hamiltonians, relevant for weak decays and mixing, are truncated at a certain 
order (typically NLO). This is not a peculiarity of the lattice approach but it depends 
on the limited number of orders which has been computed in continuum perturbation 
theory. Moreover, the choice of presenting the calculations in MS is only a traditional 
option which has not to do with the continuum or the lattice formulation of the theory. 
A discussion of this point in the case of B°-B° mixing can be found below. 

Another source of confusion is the uncontrolled propagation of values and errors of the 
parameters from one review talk to another, without verification on the original papers, 
and regardless of more recent calculations and progresses. A typical example is the value 
£ ~ 1.3. This number was only found in |38|]. All other lattice calculations find for 
this quantity values between 1.11 and 1.17 |39|-||42| in the quenched approximation and 
similar numbers have been recently confirmed by unquenched data K2| . Similarly, in . 
A. Falk quotes £ = 1.14(0.13) on the basis of a two-year old review by S. Sharpe f4"3fl . In 



the absence of unquenched calculations, the rather generous uncertainty reported in [|43f 
was justified, at the time, on the basis of theoretical estimates obtained by using quenched 
and unquenched chiral perturbation theory. These estimates have not been confirmed by 
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(partially) unquenched results, which were already available last year Irrespectively of 
recent lattice progresses, this very large uncertainty, which is not supported by any explicit 
numerical result, risks to survive and be used in future phenomenological analyses. The 



most recent figures for the S-meson decay constants, the 5-parameters, fB d yBB d and 
£ can be found in Tables |2] and [| These tables include the results presented at Lattice 



2000 JE|. In the following, for all the theoretical parameters, we have used results taken 
from lattice QCD. There are several reasons for this choice, which has been adopted also in 
previous studies of the unitarity triangle H [|, [13], [14|, [19(1 . Lattice QCD is not a model, as 
the quark model for example, and therefore physical quantities can be computed from first 
principles without arbitrary assumptions. It provides a method for predicting all physical 
quantities (decay constants, weak amplitudes, form factors) within a unique, coherent 
theoretical framework. For many quantities the statistical errors have been reduced to 
the percent level (or even less). Although most of the results are affected by systematic 
effects, the latter can be "systematically" studied and eventually corrected. All the recent 
literature on lattice calculations is indeed focused on discussions of the systematic errors 
and studies intended to reduce these sources of uncertainty. We are not aware of any 
other approach (1/N expansion, QCD Sum Rules, etc.) where such a deep investigation 
of systematic errors is being carried out for a so large set of physical quantities as in lattice 
QCD. Finally, in cases where predictions (non post-dictions) from lattice QCD have been 
compared with experiments, for example fo s , the agreement has been found very good. 

Obviously, for some quantities the uncertainty from lattice simulations is far from 
being satisfactory and further effort is needed to improve the situation. For the quantities 
considered here this is particularly true for the B®-B® mixing amplitude, as discussed 
below. Nevertheless, for the reasons mentioned before, we think that lattice results and 
uncertainties are the most reliable ones and we have used them in our study. 



4.1 Statistical and systematic effects in lattice calculation un- 
certainties 

Lattice simulations are theoretical experiments carried out by numerical integration of the 
functional integral by Monte Carlo techniques. In this respect uncertainties are evaluated 
following criteria very close to those used in experimental measurements. The results are 
obtained with "statistical errors", i.e. uncertainties originated by stochastic fluctuations, 
which may be reduced by increasing the sample of gluon-field configurations on which the 
averages are performed. It is very reasonable to assume that the statistical fluctuations 
have a Gaussian (almost Gaussian) distribution. Hence, the probabilistic inversion needed 
to infer the quantity of interest gives rise to Gaussian uncertainty models. 

To convert the results of lattice simulations in predictions for the physical amplitudes 
several steps are necessary: 

a) renormalization of the relevant operators; 

b) extrapolation to the continuum limit, namely to zero lattice spacing (a — > 0); 

We denote as partially unquenched results those obtained with two sea-quark flavours at values of 
the light-quark masses larger than the physical ones, typically of the order of the strange quark mass, or 
slightly below. 
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c) unquenched calculations. The most precise numbers have been obtained in the 
quenched approximation. Theoretical estimates and some preliminary results in 
the (partially) unquenched case are also known and they are used to estimate the 
systematic errors of the quenched results. 



a)-c) are the main sources of systematic errors for B K , fB q \j^B q an d £ and are discussed 
separately in the following. Here we want only to stress that systematic errors from 
lattice calculations are conceptually similar to some of the systematic errors present in 
experimental measurements. 

Let us consider as an example discretization errors. Since it is obviously impossible to 
work at zero lattice spacing, a method to correct for discretization effects is to compute 
a given physical quantity at different values of the lattice spacing and to extrapolate it 
to zero lattice spacing. The theory tells us whether the extrapolation has to be linear or 
quadratic in a. Thus for example one fits a given quantity Q as 

Q(a) = Q(a = 0)+Q'a + Q"a 2 + ... , (20) 

and takes Q(a — 0) as best estimate of the value of Q in the continuum limit. Since the 
"measurements" performed at fixed lattice spacing are subject to a statistical error, the 
uncertainty in the extrapolated quantity is inflated with respect to the points directly 
measured. In the quenched case, extrapolations have been made for both B^, the B- 
meson decay constants and £. Systematic studies for the 5-meson mixing parameters 
are still missing. From a comparison among calculations performed by different groups 
at different lattice spacing and with different lattice actions, an estimate of discretization 
errors can be obtained, however, also in these cases. 

When an extrapolation to the continuum limit of the lattice data has been possible, 
the final uncertainty results from the statistical error of the points measured at fixed 
lattice spacing (a residual uncertainty is present when linear or linear plus quadratic 
extrapolations give different results). Thus, in this case, it is natural to assume that the 
final error has a Gaussian distribution. 

As for the errors coming from quenched calculations, partially unquenched calcula- 
tions exist for several quantities considered in this study, namely Bx and /s ds . These 
calculations are usually performed with two light quarks in the fermion loops, at values 
of the light-quark masses larger than the physical values and an extrapolation in these 
masses is required. The calculations are generally made at a fixed value of the lattice 
spacing and thus contain discretization errors. An estimate of the quenching errors is 
obtained by comparing quenched and unquenched results at similar values of the lattice 
spacing. These comparisons are complemented by theoretical estimates of this uncertainty 



obtained by using quenched and unquenched chiral perturbation theory techniques [pHjj . 

The question arising at this point is: what is the best model to describe the theoretical 
systematic errors in lattice calculations, and hence the assessment of the uncertainty? We 
refer to the general introduction of the inferential framework given in Section |3[ As 
happens with systematic errors in experiments, this completely relies on the confidence of 
the experts about possible variations of an influence parameter, the effect of quenching 
or what would happen in passing from a perturbative order to the other or changing 
the renormalization scheme. These evaluations are unavoidably subjective, though not 
arbitrary, as long as we use the judgements of responsible experts for each input quantity. 
Using their judgements we commit ourselves too. Therefore, hereafter, when we state that 
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a parameter lies in a certain range with uniform distribution, it means that, in practice, 
we are 100% confident that the parameter lies in that region, and that, for any choice of 
a sub-interval of half the width, we are in condition of indifference (i.e. 50% confidence) 
that the value of the parameter is inside the sub- interval or somewhere else. For a more 
extended discussion see for example [45| and references therein. 



In conclusion, we cannot find any conceptual difference which would force us to treat 
experimental and theoretical uncertainties on a different footing and claim that the stan- 
dard method is a perfectly justified scientific approach able to establish confidence levels 
for the quantities of interest. In the following, for each parameter, taken from lattice QCD 
evaluations, best estimates for its central value and attached uncertainties are given. 

In order to check the stability of the results, we have also made the analysis with the 
flat part of the theoretical uncertainty increased by a factor two. 



4.2 Evaluation of the parameter fB d \jBB d 

Traditionally, the I?-parameter of the renormalized operator is defined as 

(B d \Q% B=2 (ri\B d ) = | m%jl d B Bd {^) , (21) 

where Q^ B=2 = (&7 M (1 — 75)0?) (67^(1 — 75 )d) and \i is the renormalization scale. This 
definition stems from the vacuum saturation approximation (VSA) in which B Bd = 1. 
Similarly one can define B Bs . The renormalization group invariant 5-parameter B Bd of 
Equation (H) is defined as 

B Bd = aM-rl** (l + B Bd M . (22) 

7o = 4 in all schemes whereas J depends on the scheme used for renormalizing Qf B=2 (/i). 
In the theoretical expressions, the physical amplitudes are always defined in terms of 
B Bd . The advantage is that this quantity is not only renormalization scale, but also 
renormalization-scheme independent. A residual scheme dependence remains only be- 
cause the coefficient renormalizing the lattice operator is computed at a fixed order in 
perturbation theory (NLO in this case), whereas its matrix element is computed non- 
perturbatively. This problem would arise in any approach that computes physical am- 
plitudes by combining the Wilson coefficients of the effective Hamiltonian, which are 
computed perturbatively, with hadronic matrix elements. Thus it is not specific to lattice 
calculations. The scheme dependence can be reduced by increasing the order at which 
Wilson coefficients are computed in continuum perturbation theory. 

Indeed the important quantity is not the 5-parameter itself but the combination 



(fB d yB Bd ) 2 which is used as an alias for the physical amplitude to which it is simply 
related by the factor 8m 2 B /3. This is similar to the kaon B parameter. In that case, 
however, the decay constant is taken from experiments. For the Bj meson, instead, we 
have to rely on theory also for the decay constant. Since the calculation of f Bd and B Bd 



are strongly correlated, the best way is to take the combination f Bd yB Bd from a single 
calculation rather than using f Bd and B Bd from different studies as often done in the 
literature. 
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Table 2: Decay constants from recent lattice calculations. The errors are those of the original 
publications. Some of the numbers have been taken from the recent lattice reviews [f39| [II]]. 
(*): By taking in the extrapolation the previous MILC data closer to the continuum limit only 
(corresponding to (3 = 6/<jfg > 6.0), one would find f Bd ~ 180 MeV. This is, in our opinion, 
a better extrapolation of these data. (**): This large value was found by using the Sommer 
parameter to calibrate the lattice spacing. Using the p mass, instead, they find f Bd = 186 MeV. 
We do not understand the origin of this large difference. 



The most recent calculations of the combination fB d y B Bd , in the quenching approx- 
imation, come from ref. obtained with the non-perturbatively improved action and 
a non-perturbative renormalization of the lattice operators, and from ref. ||55|| , with the 
mean-field improved action and perturbatively renormalized operators 



fB d \]B Bd = (206 ±28 ±7) MeV [51], 

fB d y/SZ = (211 ±21^7) MeV |5§. (23) 
These results correspond to the following values obtained in the same simulation 



f Bd = (174± 22+ 7 ^)MeV B Bd = 1.38 ± 0.111°;™ [51], 



<B d 

fs d 



177 ± Utll) MeV Sr, = 1.41 ± 0.061 



■0.05 
0.01 



(24) 



The numbers above agree with (our) world averages of quenched determinations, based 
on the results given in Tables |] and || which were used in |58| f\ 



fs d 



(175 ± 20) MeV 
1.36 ±0.08 ±0.05 



which can be combined to give 



B Bd = (205 ± 24 ± 8) MeV . 



(25) 



(26) 



The same quenched results for fg d is quoted in ref. 
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Table 3: i?-meson B parameters from recent lattice calculations. For comparison we have 
evolved results and errors of the original publications to a common scale. 



Note that in this case we have used decay constants and 5-parameters from different 
calculations. 

Our average of fs d does not include the results obtained with NRQCD. The reason is 
that the quenched results from two different groups are incompatible between each other: 
ref. |9l found f Bd = 147(11) (1? 2 ) (9) (6) MeV to be contrasted with the recent CP-PACS 



result |^2[ fs d = 191 (5) (11) MeV. Moreover [59| finds a rather large difference, of about 
40 MeV (from 147 MeV — > 186 MeV), between the quenched and unquenched case. This 
is not confirmed by the most recent results given in Table 0, which give differences of the 
order of 20 MeV. In our opinion, the situation with this approach is still rather confused 
and therefore we do not use the NRQCD results (until it will not be clarified). This applies 
also to the related calculations of the Bs ds parameters, which are computed within the 
same framework. 

Our average for Bs d has been computed by combining values obtained with heavy 
quark masses in the charm region, extrapolated to the B mesons (denoted as LQCD in 
Table 0), with those obtained in the HQET |57[], at lowest order in the l/m?, expansion. 
The systematic error has been estimated from the different values obtained by using only 
the LQCD results or by combining LQCD and HQET predictions. 

The world average given above for fs d includes data obtained after extrapolation to 
the continuum, as well as results obtained with improved actions at small values of a (for 
which discretization errors are expected to be smaller). A completely non-perturbative 
determination of the axial current renormalization constant Za has been also performed 
in several cases, thus eliminating this source of errors. 

In the case of the B parameter, no systematic continuum extrapolation has been 
attempted yet. In this case it is reassuring that results obtained with perturbative and 
non-perturbative renormalization techniques and for a variety of values of lattice spacing 
are so close. In the absence of any indication of large discretization errors, we ignore them 
in the following. 

As far as quenching errors are concerned, there is a general agreement that the value 
of the decay constants increases in the unquenched case. This is supported by both theo- 
retical estimates with quenched and unquenched chiral perturbation theory [Q, and by 
explicit numerical calculations, see Table |2|. The MILC and CP-PACS |54| Collabo- 
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rations find very consistent results, 



runq i rquen i i r>+0.16 

J B d I JB d — i - iz -0.11 J 

f unq /f quen = hn±QM} (27) 

respectively. Most unfortunately no unquenched determinations of the 5-parameters, or 

even better of f B \] B b , have been presented yet. For the 5-meson B parameters, chiral 
perturbation theory suggests that the unquenching error is at most of the order of 10%. 
This prediction should be supported by explicit numerical simulations which are missing 
at the moment. 

By assuming 10% quenching uncertainty for B Bd , and using the results in ( p6| ) and 
(p7|), we arrive to 

hi \fB~B d = (230 ± 25 ± 20) MeV. (28) 

which has been used in our analysis. In our preliminary analysis (mostly based on 
quenched results) which has been presented in |58| , we used instead 



f Bi ^B Bi = (220 ± 25 ± 20)MeV , (29) 
which is very close to the present value. The changes in the results for the relevant 



quantities (sin(2a), sin(2/3), 7, etc.), induced by the difference between (|28|) and (p9[), are 
negligible. 

4.3 Evaluation of the parameter £ 

Besides the B° — B° amplitude, an additional constraint is given by the ratio 



Am s \V ts \ rriB 



Am d \V td \ 2 m Bd 



r, (30) 



where ^ = f Bs \j B B J ' fB d \j B Bd - By combining the results for the decay constant ratios 
IbJ ' fs d and for the B parameters (the latter being always very close to one) (see Tables 
|2] and |3| respectively) one finds always a number of ~ 1.15. As mentioned before, there is 
no confirmation, neither in the quenched nor in the unquenched case, of a value as large 
as £ = 1.3 as found in [Q. With an estimate of the uncertainty on B B J B Bd of 10% [53 1, 
we then find 

f = 1.14 ±0.03 ±0.05 (31) 
which is the value used in our study. In our preliminary analysis we used, instead, 



f = 1.14 ± 0.06 58 . 



4.4 Evaluation of the parameter Bk 

The kaon 5-parameter, Bk-, is one of the most studied, and more accurately known, 
quantities in lattice calculations. Very precise values have been obtained, within the 
quenching approximation 

B^in = 2 GeV) = 0.63 ± 0.04 [40] 

B^(jj, = 2 GeV) = 0.62 ± 0.03 || (32) 
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which correspond to the renormalization group invariant B parameter 

B K = 0.87 ± 0.06 . (33) 

This uncertainty, meant as standard deviation, includes the contribution from the statis- 
tical error and the deviation from the extrapolation to the continuum limit (this error 
is much larger than those on individual data at fixed lattice spacing). The physical am- 
plitude has been computed at the NLO by using boosted perturbation theory [|6T|. A 



non-perturbative renormalization of the lattice operator is preferable but has not been 
performed yet. Previous experience on other quantities leads to an estimate of the error 
due to the use of the perturbative renormalization of the order of 5%. 

Uncertainties, due to the quenching approximation, have been evaluated both theo- 
retically and numerically. Using chiral perturbation theory, the error on Bk due to the 
quenched approximation has been estimated to be negligible for degenerate quark masses 
(m s = ) and of the order of 5% for realistic quark masses ||44|| . A numerical unquenched 
calculation with n/ = 2 and n/ = 4 resulted in a shift upwards of the value of Bk by 
(5 ± 2)% [|Si^l . This calculation was performed at fixed lattice spacing and only the differ- 
ence between the quenched and unquenched Bk for similar values of a was studied. Both 
the theoretical estimate and the numerical evaluation give a quenching error of (9(5%). 
The residual uncertainty is due to our ignorance on the dependence of this difference on 
the value of the lattice spacing. Since so far we have very accurate results with contin- 
uum extrapolation in the quenched case only and unquenched results without continuum 
extrapolation, it seems reasonable to consider a systematic uncertainty corresponding to 
a uniform distribution spanning the range of ±15% corresponding to the maximum dis- 
cretization effect expected in the unquenched case [fHJ. For the central value, we have 
taken the quenched result given in Equation (|33|) . We thus obtain 

B K = 0.87 ±0.06 ±0.13 (34) 

which has been used in our study. This range of values is in good agreement with others 
which can be found in the literature |63|, B4J. The allowed range for Bk in our evaluation 



and in [|63| , |64j is smaller than the one quoted in |23|, [25j, namely 0.6 < Bk < 1-0. That 
estimate tries to include results obtained with techniques different from lattice calcula- 
tions, such as QCD Sum Rules and 1/N expansion. In our opinion, the other approaches 
lack of the accuracy and control of systematic effects reached by lattice calculations for 
this parameter. We think instead that Equation (1341) corresponds to our best knowledge 



of Bk, given the present understanding of the theory. 



5 Other inputs 

In this section we briefly discuss other inputs which have been used in the present analysis: 
\Vcb\, \Vub\ an d A. These are obtained from experimental measurements combined with 
several theoretical predictions which are discussed below. 



5.1 Extraction of \V c b\ 

Using exclusive decays B® — > D* + £~ve, the value of \V c b\ is obtained by measuring the 
differential decay rate at maximum q 2 , which is the mass squared of the charged lepton- 
neutrino system. At q 2 = q^^, the D* + is produced at rest in the B^ hadron rest 
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frame and HQET can be invoked to obtain the value of the corresponding form factor 
Fd*(w = 1). The variable w is usually introduced as the product of the four- velocities of 
the B® and D* + mesons 



m 2 S o + m 2 D * 



w 



vm ■ v D *+ 

a 



2 mgo rr%D* 



w 



1 for q 2 



q 2 



(35) 



In terms of w, the differential decay rate can be written as 



dBR 



D* 



dw 



T B 487T' 



: m 3 D ,(m B - m D *) 2 K(w)Vw 2 - lF^»(w) \V C 



cb\ 



(36) 



where K(w) is a kinematic factor. 

As the decay rate is zero for w = 1, the w dependence has to be adjusted over the 
measured range. 

Four measurements obtained by the LEP collaborations have been averaged with 



the result from CLEO [p6[ , taking into account correlations induced by common sources 
of systematic uncertainties. Before averaging, results and errors have been recalibrated 
using a common set of values for the external parameters, such as the B® lifetime and 
charm hadron decay branching fractions p5|. 
The average value is 



F(l) \Vcb\ = (37.2 ± 1.5) x l(T d , % l NDF = 9-15/4 



(37) 



The "fit probability"^] is only 6%, giving rise to the suspect that there could be some extra 
systematic effect playing an important role. 

Another evaluation of this average has thus been done following a model developed 
to combine results which appear to be in mutual disagreement. It consists in assuming 
that quoted uncertainties (sj) for each measurement are proportional to the unknown real 
uncertainty (o"j/rj = Si). A distribution probability for is then assumed. A simple 
model, which depends on two parameters, 5 and ( is given by 



f(r) 



2 ^ r -(25+i) eX p-CA 



(3* 



The natural choice for S and ( corresponds to an expected value of 1 for r^, with 100% 
uncertainty. The final results are largely independent from the precise value of the pa- 
rameters. We have chosen the values of ( = 0.6 and 5 = 1.3 taken from Ref. |69 . 



The probability distribution for the quantity F(l) \V c b\ is then obtained using the 
Bayes theorem 



f(x) oc 



exp 



1 5 

9 H ( x ~ Xi)w itj (ai/ri, aj/rj)(x - Xj) 



Y[f(n)dri 



(39) 



\i=l 



In this expression, x = F(l) \V c t,\ and Wij is the weight matrix. The latter is obtained by 
inverting the error matrix, after the inclusion of the uncertainties given by the quantities 

6 "p-value" would be the correct modern statistics term to be used instead of "fit probability" or "x 2 
probability" p7i 68j. These expressions can be highly misleading as explained in |p4|. 
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Ti (considered to be independent for the different measurements). From this distribution, 
which has non-Gaussian tails, the average and the standard deviation have been obtained 

F(l) \V cb \ = (37.1 ± 1.9) x 1(T 3 . (40) 

The central value is practically the same as obtained with the usual fit and the standard 
deviation has increased by 30%. 

From the exclusive measurements and using F(l) = 0.88 ± 0.05 J7[J, \V c b\ has been 
obtained 

\V cb \ = (42.2 ±2.1 ±2.4) x 10~ 3 . (41) 



-1 




Vcb*10 3 



Figure 1: Probability distribution function for \V c b\ obtained by combining exclusive and 
inclusive measurements of 6-hadron semileptonic decays. The full line distribution corresponds 
to a Gaussian whereas dotted curves are obtained using the procedure explained in the text. 
Two sets of values for the £ and 5 parameters, namely (0.6, 1.3) and (1.4, 2.1) have been 
tried and correspond to almost identical distributions (dotted curves). The quoted error bars 
for individual results correspond to all sources of uncertainties but the combination procedure 
accounts for correlations, coming from common systematic uncertainties. 

The inclusive measurements of the semileptonic branching fraction of 6-hadrons give 
instead 

|y c6 | = (40.7 ±0.5 ±2.0) x 10~ 3 . (42) 

The interested reader may consult |70| for more details on the quoted uncertainty for 
inclusive decays. 
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Figure 2: The p.d.f. for \V ub \. 



The combination of all results, using the procedure explained previously when aver- 
aging F(l) \Vcb\ measurements, taking into account correlated systematics, gives 

(41.0 ± 1.6) x icr 3 . 



\v. 



cb\ 



(43) 



The corresponding p.d.f, which has been used in the present analysis, is shown in Figure 



5.2 Extraction of \V u b\ 

\V u b\ has been obtained from the measurements done at CLEO and LEP. The CLEO 
collaboration |71] has measured the branching fraction for the decay B% — > p + £~u e and 
deduced a value for \V u b\ using several models to describe the decay form factors. LEP 



collaborations JF2] have developed dedicated algorithms to be sensitive to a large fraction 
of the inclusive decay rate b — ► u£~v e and, with some assumptions, a value for \V u b\ is 
obtained. The two measurements are 



(32.5 ± 2.9 ± 5.5) x 10" 4 CLEO 
V uh \ = (41.3 ± 6.3 ± 3.1) x 10" 4 LEP 



IK&I = 



(44) 



where the second uncertainty is theoretical. The p.d.f. for the CLEO measurement 
is thus a convolution of a Gaussian and a flat distribution. The theoretical error for 
the LEP measurement, being the convolution of several different errors, is taken from a 
Gaussian distribution [72] . Combining the two distributions, in practice we obtain almost 
a Gaussian p.d.f. (see Figure H) corresponding to 



|K&| = (35.5 ±3.6) x 10" 



(45) 
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5.3 The parameter A 
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Figure 3: Probability distribution function for A obtained by combining measurements of \ V u d\ 
and \V US \. The full line distribution corresponds to a Gaussian whereas the dotted curve is 
obtained using the procedure explained in the text. 



Measurements of \V u d\ and |V^ S | reported in [3lJ have been combined f\ using the pro- 
cedure explained in Section [54], assuming that \V u d\ = cos 6^ and |V^ S | = sin# c = A. The 



additional contribution from \V u b\ to the unitarity condition can be safely neglected owing 
to present uncertainties. The p.d.f. for A strongly deviates from a Gaussian (see Figure 
H) because present determinations of this quantity from \V uc t\ and \V US \ measurements are 
more than two standard deviations away. The average and the standard deviation of the 
obtained distribution correspond to 

A = 0.2237 ±0.0033 (46) 

whereas a classical fit gives A = 0.2225 ± 0.0019. 

Using the values of \V c b\ and A given respectively in Equations (f43|) and (^) the following 
result is obtained for the parameter A 



.4 



W, 



eft 



A 2 



0.819 ±0.040 



(47) 



It may be noticed that, in the present determination of the p and r\ parameters, \V c b\ and 
A, but not A, are used as constraints. 

7 The corresponding values used in the present analysis are : \V u d\ = (0.9740 ± 0.0010) from nuclear 
beta decays, \V ud \ = (0.9728 ± 0.0012) from neutron decay and \V US \ = (0.2196 ± 0.023) from K e3 decays. 
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6 A new procedure to include the information com- 
ing; from B° — B° oscillation analyses 

° s s J 

In this Section, we illustrate a new procedure to include in the analysis the experimental 
searches for B° s — mixing. It makes use of the likelihood along the lines discussed in 
Section [3]. 

The study of B° s — B® oscillations is done by introducing the oscillation amplitude A as 
explained in [73]. Events have been distributed in two classes corresponding to oscillating 
and non-oscillating candidates, according to the information obtained by tagging the 
presence of a B or a B meson at the beam interaction time and at the 6-hadron decay 
time. The expected decay-time distribution of events is obtained by using a function 
which contains time distributions for all components, taking into account the non-perfect 
event classification, the different behaviour of 6-hadron decays, which depend on their 
type {B®, B~ , B® or 6-baryon), and background components originating from charm and 
light flavours. The time distribution for the oscillating (non-oscillating) signal is given by 



V(B» -> B»(B»)) 



2t b o 



[1 - (+)AcosAm s t)e 



(4f 



Such theoretical time distributions have been convoluted with the expected time reso- 
lution of the measurements. For each value of Am s and each analysis, the oscillation 
amplitude A is fitted. Fitted amplitude values have been combined using a common 
set of determinations for external parameters and taking into account possible correla- 
tions between the different analyses. In this respect the use of the oscillation amplitude 
provides a simple framework to account for all these effects. The 95% CL. limit on Am s 
corresponds to the value Am^ m for which A + 1.645 er^ = 1, whereas the sensitivity 
Am s s ens corresponds to 1.645 o A = 1. It is the value of Am s at which it is expected that 
the 95% limit will be set, if A = 0. In addition, the final likelihood distribution can be 
retrieved from the combined amplitudes corresponding to all analyses as explained in the 
following. 

The information coming from oscillation amplitude measurements, which gives con- 
straints on possible values of Am s , was included [|1|, [b|, [15], [23], ^4j, up to now, using the 
distribution PI 

1 /A- 1\ 2 

F(Am s ) = exp--(- — -) . (49) 



2 V <r A J 

Recently it has been proposed to use the log-likelihood function Alog£ 



erenced to its value obtained for Am s = oo 
a different context, have been detailed in 
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(Am s ) ref- 
|. Similar considerations, developed in 
The log-likelihood values can be easily 



deduced from A and a A using the expressions given in [73 



Alog£°°(Am s 



A-l 



Alog£°°(Am 
Alog£°°(Am 



s ) nomix 



A 



1 

2 ~ 
1 1 



2aY 



1 1 

2a\ 



The dependence of A and a a on Am s is always implicitly assumed. 



a 



2 ' 
A 



(50) 
(51) 
(52) 



21 



The last two equations give the average log-likelihood value when Am s corresponds to 
the true oscillation frequency (mixing case) and when Am s is far from the oscillation 
frequency (\Am s — Am* rue | ^> T/2, no-mixing case). T is the full width at half maximum 
of the amplitude distribution in case of a signal; typically T ~ l/r B o. In the following it 
is shown that: 

• Equation (ff9|) does not represent the optimal way to include the Am s information. 
In particular it is incorrect in case of a measurement; 

• Equation fl5Ti| ) allows to define the likelihood ratio R 

R(Am s ) = e -Alog£~(Am a ) = C(Am a ) 

C(Am s = oo) 

The function R corresponds to the ratio of probability densities for different Am s 
values. The absolute Am s probability density, instead, remains undefined because 
C(Am s ) stays constant for Am s values much larger than the sensitivity; in this 
region R is equal to unity. 

• a reasonable procedure to extrapolate R in regions where the direct measurements 
of the amplitudes are not available can be established. 

6.1 Comparison between the new and the old methods 

The main concern, in the previous approach, was that the sign of the deviation with 
respect to the value A = 1 was not used, whereas it is expected that an evidence for a signal 
would manifest itself by giving an amplitude value which is simultaneously compatible 
with A = 1 and incompatible with ^4 = 0. 

The new method, at variance, includes the relative weight of the two hypotheses. 
Equations (|50|) and ([k|) show that R has a clear interpretation as ratio between two 
cases: mixing (A = 1) or no mixing (^4 = 0). The properties of the log-likelihood ratio, 



R, and its application in Higgs boson searches have been introduced in [75 . 
Another problem of the previous approach is that the sign of the deviation of the amplitude 
with respect to unity is not considered: this implies that a lower probability is attributed 
to Am s values with A > 1 with respect to Am s values having A = 1. Since this behaviour 
is clearly undesired, in the old method the amplitudes larger than unity have been set to 
unity (as it has been implicitly done in Hl3|-p3|). 



The difference between the old and the new method are illustrated in two cases: 

• a single analysis having the sensitivity at 15 ps _1 (see next paragraph and Figure 
H) simulating a signal at Am s = 5ps _1 (this value has been chosen well below the 
sensitivity in order to amplify the effect). 

• the world average analysis (see Figure |5]) 

In Figure [§-a and Figure |]-b, respectively, the amplitude spectrum A and the corre- 
sponding distribution for A£°°(Am s ) are given for the simulated case, in which there is 
a signal at Am s = 5 ps _1 . Figure ^-c shows that there is a marked difference between 
the behaviour of R and F. While the function R shows, by construction, a maximum 
corresponding to the fitted value (Am s = (5.10 ± 0.25) ps _1 ), the function F is not able 
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Figure 4: Amplitude analysis with Am s generated at 5 ps _1 : a) amplitude spectrum, b) 
A log £°°(Am s ), c) comparison between R and the function F used in the previous method. 
The maxima of the two functions have been rescaled to unity in order to ease the comparison 
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Figure 5: World average amplitude analysis: a) amplitude spectrum, b) Alog£°°(Am s ), c) 
comparison between R and the function used in the previous method. 

to spot the maximum, in particular it attributes large probability to any Am s value with 
A ~ 1 regardless its error. 

For the world average analysis, Figure |5] shows that the agreement between R and F is 
acceptable only when A ~ 0. 
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6.2 How to treat Am s regions without amplitude measure- 
ments? 

The procedure used to continue R beyond the last measured amplitude value (Am' ost ) is 
based on the continuation of the amplitude spectrum from which R is deduced: 

• continuation of a (A): the behaviour of cr(A) for Am s > Am|, ,m can be reproduced 
by tuning the parameters of a fast simulation (toy-MC). The method used here is 



similar to the one presented in |7q| . The errors on the amplitude can be written as: 
a -\A) = VNr]B a (2e - 1) W(a L , a P , Am s ) 

where N is the total number of events, t]b s the purity of the sample in B s decays, e 
the purity of the tagging at the decay time, ap is the B s flight length uncertainty 
and ap the relative uncertainty on its momentum. The parameters <jl, <Jp and the 
global factor that multiply W have been obtained by adjusting the simulated error 
distribution on the one measured with real events. 

Figure |6| shows the agreement between the toy-MC calculation and real data for 
Am s < Am'" 1 (the upper bound on Am s of the amplitude plot for which there are 
measurements) and the extrapolation to higher values. It has been verified (with an 
independent toy-MC) that the extrapolated errors approximate well the simulated 
error distribution. 

• continuation of A: this part is more critical than the previous one. In particular it is 
more sensitive to the real amplitude spectrum. Nevertheless if Am s s ens << Ara^ 1 , 
the significance S (S = A/aX) is approximately constant. 

In the following the amplitudes at Am s > 25 ps _1 (which corresponds to Amf for the 
combined world average plot) are deduced using the extrapolation of the amplitude errors 
and fixing the significance at the last measured value. As a consequence, for Am s values 
larger than Am' asi , R is fixed to the last measured point. Although this procedure is 
reasonable, it should be stressed that it is very desirable to have all the amplitudes (with 
errors) up to the Am s value where R approaches its plateau. 



6.3 The likelihood function R used in the analysis 

The present world average on Am s shows an "hint" of signal at Am s ~ 17.5 ps" 1 (see 
Figure [5|-b). Its significance (about 2.2 a) is not sufficient to claim evidence for an 
observation of the B° s — oscillation. This effect is in fact larger than the average value, 
of ~ 1.2 cr, which is expected in case of a real signal in that region. 
Whatever conclusion can be drawn from these data, it should be stressed that: 

• the likelihood function R deduced from data contains the present knowledge on Am s 
and it is the optimal weight function to be used in the fit; 

• the use of R does not imply any assumption either on the evidence or on the presence 
of a signal. It only translates the experimental fact that the log-likelihood function 
has a minimum at 17.5 ps -1 with a 2.2 o significance. 
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Figure 6: Comparison between the error distribution computed with a toy-MC (solid line) and 
the measured amplitude errors (circles). 




Figure 7: Comparison between different weight functions used in the fit. The curve determined 
by the continuous and dashed lines corresponds to the likelihood ratio obtained with, respec- 
tively, the measured amplitudes and the extrapolation to large values of Am s . The dotted curve 
is obtained by truncating the previous distribution once the amplitude becomes larger than one. 
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Although we believe that the use of R is, indeed, the best method to include the 
information on Am s , we briefly mention some alternative approaches. The main difference 
is that, with these methods, a fraction of the available experimental information is lost. 

1. 95% CL. limit: this method uses a stepwise function starting at Am^ m . We mention 
it only for completeness since the choice of the CL. is arbitrary. 

2. Min(R,l): in this approach it is considered that no value of Am s has to be preferred 
with respect to Am s = oo (only exclusion is possible). It was used in our prelimi- 



nary analyses |58[ and consists in setting A = 0.5 for all amplitudes having larger 
values (see equation (|50|) . Since there is no reason to throw away the information 
contained in amplitudes with values ranging between 0.5 and 1, this method has 
been abandoned. 

3. Min[R(v4 = 1),R]: in this method one tries to avoid biases induced by "lucky" 
fluctuations. The strategy is then to limit the likelihood ratio to the value obtained 
when the amplitude reaches A = 1 for the first time. The effect is to flatten the 
ratio R around its maximum. 

Two different definitions of R have then been used as weight function in the fitting 
procedure: 

• the complete R distribution deduced from data and from the toy-MC; 

• the function Min[R(A = 1),R]. 

They are shown in Figure |^ and the second choice has been considered as a "pessimistic" 
approach to evaluate the induced variation on the fitted quantities in Section |[ 



7 Results and discussion 

In this Section we give the results for the quantities of interest. Values of the hadronic 
parameters, which can be determined quite accurately by assuming the validity of the 
Standard Model, are also extracted. The central value is always given using the average, 
and the error corresponds to the standard deviation. For asymmetric p.d.f. we also give 
the median and the error corresponds to regions containing 34% probability on each side 
of the median. 

7.1 Results obtained with all measurements 

The region in the (p, ff) plane selected by the measurements of \ek\, |V^/V^,|, Am d and 
from the information on Am s (using the R function of Figure [7|) is given in Figure |8|. In 
Figure || the uncertainty bands for the quantities, obtained using Equations dH)-(§), are 
presented. Each band, corresponding to only one of the constraints, contains 68% and 95% 
of the events obtained by varying the input parameters. This comparison illustrates the 
consistency of the different constraints provided by the Standard Model. In the present 
studies, two statistically equivalent procedures have been used. They differ in the way 
values for p and f] have been extracted and in the inclusion of QCD corrections to ex and 
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Figure 8: Allowed regions for p and fj using the parameters listed in Table 0. The contours at 
68%, 95%, 99% and 99.9% probability are shown. The full lines correspond to the central values 
of the constraints given by the measurements of \V u b\ / \Vcb\, \^k\ and Amj. The two dotted 
curves correspond, from left to right respectively, to the 95% upper limit and to the value of the 
sensitivity obtained from the experimental study of B® — B® oscillations. 




P 

Figure 9: The allowed regions for p and f\ (contours at 68%, 95%) are compared with the 
uncertainty bands (at 68% and 95% probabilities) for \V U \\ j\y c b\-, \^k\, Am^ and the limit on 
Am s /Amd (dotted curve). 
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to B° — B° mixing, which are either computed [JIS or taken as independent inputs |TJ . 
The measured values of the two parameters are 



and 



p = 0.224 ±0.038, 



0.221 ±0.037, 



f) = 0.317 ±0.040 



0.315 ±0.039 



(54) 



(55) 



using the methods of ref. |14]] and [19| respectively. The two quantities are practically 
uncorrelated (correlation coefficient of -5%), as it can be seen from the contour plot of 
Figure |8|. Fitted values for the angles of the unitarity triangle have been obtained also 

sin(2/3) = 0.698 ±0.066, sin(2 a) = -0.42 ± 0.23 , 7 = (54.8 ± 6.2)° [14] (56) 

(the quoted error is symmetric, in spite of the small asymmetry of the distributions shown 
in Figure [U]) and 

sin(2/3) = 0.692 ± 0.065 , sin(2a) = -0.43 ± 0.21 , 7 = (54.9 ± 5.7)° [19]. (57) 
respectively. 

The results in the two approaches are in agreement and, in the following, quoted values 



and figures are those obtained using the approach of [Q. In Figure [K| the p.d.f. for the 
angles of the unitarity triangle are given. Few comments can be made: 

• sin(2/3) is determined quite accurately. This value has to be compared with re- 
cent measurements of this quantity using J/ipK s events from LEP [fT^l , CDF |78fl , 
BaBar |7j| and BELLE [5D|. The accuracies of these measurements are completely 
dominated by the statistical errors and their average is sin(2/3) = 0.49 ± 0.16. 

• the angle 7 is known within an accuracy of about 10%. It has to be stressed that, 
with present measurements, the probability that 7 is greater than 90° is only 0.03%. 
Without including the information from Am s , it is found that 7 has 4% probability 
to be larger than 90°. The central value for the angle 7 is much smaller than that 
obtained in recent fits of rare I?-meson two-body decays |H|]. It remains to be seen 
to which extent the results of [|5T| are affected by the model dependence in the 
theoretical description of two-body decays. In this respect, it would be interesting 
to examine under which conditions these decays can be described by the same value 
of 7 as found in the present study. An exploratory work in this direction can be 



found in 82 



7.2 The CKM triangle from 6-physics alone 



As four constraints are used to determine the values of two parameters, it is possible to 
relax, in turn, one (or more) of these constraints, still obtaining significant confidence 
intervals. An interesting exercise consists in removing the theoretical constraint for Bk 

|). The corresponding selected region in the (p, ff) 



in the measurement of \sk\ (|18 



plane is shown in Figure [11], where the region selected by the measurement of \£k\ alone 
is also drawn. This comparison shows that the Standard Model picture of CP violation in 
the K system and of B decays and oscillations are consistent. In the same figure, we also 
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Figure 10: The p.d.f. for sin(2a), sin(2/3) and 7 . The red (darker) and the yellow (clearer) 
zones correspond respectively to 68% and 95% of the normalised area. 
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Figure 11: The allowed regions (at 68%, 95%, 99% and 99.9% probability) for p and fj using 
the constraints given by the measurements of \ V u b \ / \ Vcb\, Am^ and Am s . The constraint due to 
\ek\ is not included. The regions (at 68% and 95% probability) selected by the measurements of 
\ek\ (continuous (la) and dotted (2a) curves) and sin(2/3) (darker (la) and clearer (2a) zones) 
are shown. For sin(2 (3) the two solutions are displayed. 



compare the allowed regions in the (p, fj) plane with those selected by the measurement 
of sin(2/J) using J/ipK s events. 

Using constraints from 6-physics alone the following results are obtained 

fj = 0.302+°:oe? ! [0.145 - 0.400] at 95% ; > 0.08 at 99% 
sin(2 (3) = 0.678lg;5o? ; [0.392 - 0.818] at 95% (58) 

(in terms of average and standard deviation the results are fj = 0.296 ± 0.063 and sin(2 /3) 
= 0.663 ± 0.109). 

Another way for illustrating the agreement between K and B measurements consists in 
comparing the values of the Bk parameter obtained in lattice QCD calculations with the 
value extracted from Equation (§), using the values of p and fj selected by 6-physics alone 

^b-phys = 0.90+0.30^ Q g4 ^ ^b-phys < Lg ^ ^ % probabiUty (59) 

(in terms of average and standard deviation the result is _B^ _phys = 0.97 ± 0.23). 

Since Bk is not limited from above, for the present study, probabilities are normalised 

assuming Bk <5. The p.d.f of B K ~ phys is shown in Figure [12]. 

It can be noticed that the values of two theoretical parameters, namely fs d \J Bs d and 
£ have been used to obtain this result. Hopefully, when accurate direct measurements of 
sin(2 (3) will become available, it will be possible to remove another theoretical input (or 
even all of them). 
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Figure 12: The p.d.f. for B a K 



b— phys 



7.3 The expected value for Am s 

Figure shows the allowed region for p and 77 obtained when removing the constraint 
coming from the study of B° s -B° s mixing. It illustrates the importance of the use of this 
information. This can be illustrated also, from the p.d.f. of the angle 7 obtained with or 
without including the Am, constraint (see Figure [H]). High values for 7 are excluded at 
high confidence level by the experimental lower limit on Am,. 

It is also possible to extract the probability distribution for Am s , which is shown in 
Figures |15| (see also |8l|). From this distribution one obtains 



Am, = (16.3 ±3.4) ps 



-1 



9.7 < Am, < 23.2 ps" 1 at 95% probability. (60) 



If the information from the B° s — B° s analyses is included, results become 

Am s = (17.3+^) ps" 1 15.6 < Am, < 20.5 ps" 1 at 95% probability 



(61) 



(in terms of average and standard deviation the result is Am s = (17.7 ± 1.3)ps _1 ). 
These values are in agreement with the recent estimate of Am, = 15.8(2.3)(3.3) ps" 1 , 
presented in |5l|. 



7.4 Determination of fB d yBB 



The value of /s d y-B_B d can be obtained by removing the theoretical constraint coming 
from this parameter in B®-B® oscillations. Using the two other theoretical inputs, Bk 
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p 

Figure 13: The allowed regions for p and fj using the constraints given by the measurements 
of \£k\, \V u b\ I |K& I an d Am^ at 68% and 95% probability are shown by the thin contour lines. 
Selected regions for p and rj when the constraint due to Am s is included have been superimposed 
using thick lines. The different continuous circles correspond to fixed values of Am s . Dashed 
circles, drawn on each side of the curve corresponding to Am s = 18.0 indicate the effect 

of a variation by ±0.06 on £. 



and £ , fB d \jB Bd is measured with an accuracy which is better than the current evaluation 
from lattice QCD, given in Section 0. From the p.d.f. shown in Figure 16L 



f Bd ^B Bd = {2?>l±lh) MeV 



(62) 



is obtained. The present analysis shows that these results are in practice very weakly 



dependent on the exact value taken for the uncertainty on f Bd yB Bd . An evaluation 
of this effect has been already presented in |15| where the flat part of the theoretical 



uncertainties on f Bd \l B Bd was multiplied by two. Similar tests will be shown in Section 



7.5 Further comments on the Lattice predictions for fB d \/BB d 
and Bk 



The values found in the present analysis for f Bd y B Bd and Bk are in agreement with those 
of previous studies [Il|-[50|. They are also in agreement with the predictions from lattice 
QCD. It can be noticed that lattice predictions for these quantities existed well before 
it were possible to extract them from the analysis of the unitarity triangle. Moreover 
these predictions have been stable over the years: one of the first calculations of Bk gave 
5]f = 2 GeV) = 0.65 ±0.15 [HJ in 1987 (corresponding to B K = 0.90 ± 0.20) and B K 
was estimated to be Bk = 0.84 ± 0.03 ± 0.14 in 1996 j3; a compilation by one of the 
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Figure 14: The p.d.f. for 7 obtained with (red-clearer) and without (blue-darker) including 
the Am s constraint. The vertical line corresponds to 7 = 90°. The probabilities of 7 > 90° are 
also given. 
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Figure 15: Probability distribution of Am s . The information from B® — B® oscillations is 
used(not used) (distribution on the right(left). 
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Figure 16: The p.d.f. for f Bd jB Bd 



authors of the present paper gave f Bd 
MeV, in 1995 and 1996 respectively fl§6 



B B d 



(220 ±40) MeV and f Bc 



B B d 



(207 ±30) 



7.6 Lower bounds on fB d \l^B d and Bk 



The region in the plane (/b^v-Bb , Bk), which is obtained by removing the theoretical 



constraints on these quantities, is shown in Figure \L7\. It appears that present constraints 
can cope with very large values of Bk and thus it was needed to restrict the possible range 
of variation for this parameter. In the present study, probabilities have been normalised 
assuming Bk < 5. The 68% and 95% probability contours have rather different shapes. 

Within 68% probability, both f Bd J B Bd and Bk are well constrained. The most important 
conclusion which can be drawn from this study is the simultaneous lower bounds on 



f Bd \JB Bd and B K , namely 



B K > 0.5 and /„ 



B Bh > 150 MeV . 



at 95% probability. 



(63) 



7.7 Other quantities of interest 



For completeness the values of the phases x, x' an d Im X t |[7 are also given 
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Figure 17: The 68% and 95% contours in the (fB d y Bb^Bk) plane. 



v* d v us \ ./ (l^l 2 ^) 



X = arg =atan 



V c * d V cs J \l- A 2 /2 - (\V cb \ 2 /2)(1 - 2p) 



ImA t = A \V cb \ 2 fj. (64) 
The measured values of these parameters are 

X = (17.4 ±2.1)°, x' = (0-031 ±0.003)°, Im X t = (1.19 ± 0.12) 10~ 4 . (65) 



8 Stability of the results 

The sensitivity of present results on the assumed probability distributions attached to 
the input parameters was studied. The comparison of the results obtained by varying 
the size of the theoretical uncertainties has been done to evaluate the sensitivity to these 
variations of uncertainties quoted on fitted values. This must not be taken as a proposal 
to inflate the uncertainties obtained in the present analysis. 

In these tests, all values for uncertainties of theoretical origin have been, in turn, 
multiplied by two. For the quantities \V u b\ and \V c b\, new p.d.f. have been determined, 
following the prescriptions given in Sections |5.1| and |5.2| and used in the analysis; central 
values and uncertainties quoted in Tables [| and [| correspond to the average and standard 
deviation of these distributions. 

The main conclusion of this exercise is that, even in the case where all theoretical 
uncertainties are doubled, the unitarity triangle parameters are determined with an un- 
certainty which increases only by about 1.5 (see Tables ^ and ^). 
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Parameter (s) 


Modified Value(s) 


p 


V 


All 


None 


0.224 ±0.038 


0.317 ±0.040 


B K 


0.87 ±0.06 ±0.26 


0.225 ±0.040 


0.314 ±0.045 


fB d y B Bd 


(230 ± 25 ± 40)MeV 


0.224 ±0.039 


0.318 ±0.041 


t 


1.14 ± 0.04 ± 0.10 


0.221 ±0.045 


0.318 ±0.042 


\V ub \ 


(37.8 ± 4.9) 10~ 4 


0.240 ±0.044 


0.334 ±0.048 


\Vcb\ 


(40.7 ± 2.6) 10- 3 


0.224 ±0.039 


0.320 ± 0.048 


Am s cons. 


see Sect. || 


0.227 ±0.047 


0.315 ±0.042 


All 


as above 


0.249 ±0.064 


0.333 ±0.065 



Table 4: Stability tests (1). Central values and uncertainties for p and n have been defined 
relative to the median of the corresponding p.d.f. 



Parameter (s) 


Modified value(s) 


sin (2/3) 


sin (2a) 


7 (degrees) 


All 


None 


0.698 ±0.066 


-0.42 ±0.23 


54.8 ±6.2 


B K 


0.87 ±0.06 ±0.26 


0.693 ±0.071 


-0.45 ±0.26 


54.5 ±6.7 


/s d \J B B d 


(230 ± 25 ± 40)MeV 


0.698 ±0.067 


-0.44 ±0.24 


54.9 ±6.3 




1.14 ± 0.04 ± 0.10 


0.697 ±0.067 


-0.42 ±0.26 


55.4 ±7.2 


\v ub \ 


(37.8 ± 4.9) 10~ 4 


0.733 ±0.080 


-0.40 ±0.24 


54.4 ±6.1 


\v cb \ 


(40.7 ± 2.6) lO- 3 


0.701 ±0.077 


-0.43 ±0.26 


54.8 ±6.6 


Am s cons. 


see Sect. ^| 


0.698 ±0.067 


-0.45 ±0.27 


54.6 ±7.6 


All 


as above 


0.738±g:SS 


-0.42 ±0.37 


53.3 ± 10.0 



Table 5: Stability tests (2). Central values and uncertainties for sin (2/3), sin (2a) and 7 have 
been defined relative to the median of the corresponding p.d.f. 



9 Comparison between the standard and the 
95%C.L. scanning approaches 

The theoretical basis which allows to define, in the standard approach, regions of the 
(p, 77) plane which correspond to any given value of confidence has been already discussed. 
For completeness, in this Section, a comparison of present results at the 95% CL. with 
the corresponding regions selected using the 95% CL. scanning approach is made. In the 
latter approach, this region is the envelope of 95% CL. contours denned at several points 
and, as discussed previously, it is not easy to understand to which level of confidence they 
correspond to. 



The results of the study given in |24j have been used in this comparison and the 
same central values for the parameters have been used in the two cases. When, in the 
95% CL. scanning approach, a parameter is scanned over a given interval, in the standard 
approach a flat probability distribution defined over the same range has been used. In 
Table ^| central values and uncertainties used in this comparison, taken from [23], have 
been collected. 

The constraint coming from the study of — B° s oscillations is included using the 
weight given in Equation fl49|) , since it was used also in the 95% CL. scanning analysis. 

In Figure the "95% CL." contours obtained with the two methods have been 
compared and the 95% CL. intervals for p and 77 are given in Table [7]. In Figure |T9|-left 
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parameter 


Value ± Gaussian ± Flat errors 


\Vcb\ 


(40.0 ±2.0 ±0.0) x 10~ 3 


\Vub\ 


(34.0 ± 1.4 ±5.0) x 10" 4 


i — * ' > L CL 


(0.473 ± 0.016 ± 0.000) ps" 1 


Am s 


> 14.3 ps" 1 


\ £ k\ 


(2.285 ± 0.018 ± 0.000) x 10" 3 


B k 


0.80 ±0.00 ±0.15 


fB d \jB Bd 


(200 ± ± 40) MeV 


C _ /flj V ^B s 


1.14 ±0.00 ±0.08 







Table 6: Values of the quantities used in the comparison between the standard and the 
95% CL. scanning methods. The quoted flat error values are used to define the range 
of variation of the corresponding parameter in the 95% CL. scanning approach and to 
define the uniform probability distribution for this parameter in the standard method. 



parameter 


standard approach 
(95% prob. range) 


95% CL. scanning approach 


P 
fj 

sin(2 0) 


-0.06, 0.31 
0.26, 0.42 
0.56, 0.82 


0.00, 0.30 
0.20, 0.45 
0.50, 0.85 



Table 7: Comparison of standard and 95% CL. scanning approaches. Note that intervals 
quoted for the results of the standard approach correspond to 95% probability for each quantity, 
whereas in the 95% CL. scanning method the intervals are obtained from the maximum and 
minimum values spanned by the 2-D envelope. 



the distribution of p is also shown. 
Several remarks are in order: 

• when using the same values for the input parameters and for their uncertainties, the 
95% CL. scanning and the standard approaches select similar regions in the (p, fj) 
plane. It is then unjustified to qualify as too optimistic the standard approach. The 
main difference is precisely the opposite, namely that, unlike the case of 95% CL. 
scanning, with the standard method one is able to quantify the confidence corre- 
sponding to any given interval. 

• in the standard approach the probability distribution in the (p, fj) plane is mean- 
ingful. In the present exercise, the interval [0.15, 0.30] for p is more probable than 
[0., 0.15]. Quoting simply that the value of p is expected to be between and 0.3, 
at 95% CL. does not indicate that there is a most probable value around 0.2. Such 
considerations are not in order in the 95% CL. scanning approach; 

• even in the 95% CL. scanning approach it is necessary to define intervals which are 
expected to contain the true values of the parameters. But then some probability 
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Figure 18: The contour corresponding to 95% probability in the standard approach, given 
by the external continuous line, has been compared with the region spanned by 95% C.L. 
ellipses obtained in the 95% C.L. scanning method. The internal contour, corresponding to 68% 
probability in the standard approach, is also drawn. 
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Figure 19: Left: Probability distribution of the p variable obtained in the standard approach. 
The interval selected in the '95% C.L. scanning method has been indicated. Right: Compari- 
son between different input distributions for the parameter fs d \J L>B d an d the final probability 
distribution obtained with all constraints but the measurement of Am^. 
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has to be attributed to these intervals. It is thus illusory the statement that with 
this method one does not attribute any probability distribution to the theoretical 
inputs. Indeed the statement that an input cannot assume a value outside a given 
interval, is equivalent to attribute zero probability to the region outside the interval; 

the flat distribution for a theoretical parameter, in some cases, is a very strong 



a-priori, as illustrated in Figure [19|-right, where the choice of the central value 

and error on fs^Bs used in [2~1| is shown. The flat interval corresponds to the 
values quoted in Table The broad Gaussian-like distribution, used in the present 
analysis, results from a convolution between a flat and a Gaussian distributions 
using the parameters of Table |2|. This is also shown in the figure, together with the 
final p.d.f. resulting from our analysis. From this comparison, it appears that the 
region favoured by present measurements is marginally compatible with the interval 
selected in the 95% C.L. scanning. In particular, the latter excludes a priori (large) 

values of fs\l Bb which are still compatible with data. 



10 Conclusions 

In this paper, we have discussed the inferential framework which allows a consistent 
analysis of the unitarity triangle, where both experimental and theoretical uncertainties 
play an important role. In particular, we were concerned with the treatment of the 
theoretical errors, which raised several controversies in the past. We have shown that, 
in this framework, there is a consistent way of handling the theoretical uncertainties and 
that these can be included in the analysis in a way similar to the experimental ones. The 
main results of our analysis have been discussed in Section [7] and, for convenience, are 
listed below (see Table ||). 

Present measurements of \sk\, \V u b/V c b\, Arrid and of the limit on Am s allow, within 
the framework of the Standard Model, and after including results from lattice QCD on 

Br, fB d \J Bs d and £, to determine the parameters of the CKM unitarity triangle. 



Parameters 


Present analysis 


P 


0.224 ± 0.038 


fj 


0.317 ± 0.040 


sin(2 0) 


0.698 ± 0.066 


sin(2a) 


-0.42 ± 0.23 


7 


(54.8 ± 6.2)° 



Table 8: The results on the unitarity triangle parameters presented in this paper. 



The dependence of central values and quoted uncertainties on values and uncertainties 
assumed for the theoretical parameters, and on different procedures to include the avail- 
able information, have been evaluated. These exercises demonstrate that the standard 
approach employed to get our numbers is on a sound basis and that results and uncer- 
tainties for the quantities of interest are stable. The determination of central values and 
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uncertainty distributions for the theoretical inputs used in the present analysis has been 
reviewed and explained. 

By removing, in turn, the different constraints it has been verified that, within the 
present uncertainties, the Standard Model description of CP violation is in agreement 
with the measurements. In particular, B oscillations and decays select, in the (p, ff) 
plane, a similar region as \sk\ f° r the K° system. This constitutes already a test of com- 
patibility between the measurements of the sides and of the angles of the CKM triangle. 
Its accuracy is given by the quoted uncertainties on Bk appearing in the first row of Table 
1 

This Table gives a summary of the determination of the values of some theoretical pa- 
rameters and of Am s obtained in the present analysis and corresponding estimates from 
lattice QCD: 



Parameters 


Present analysis 
(with the contraint) 


Present analysis 
(without the contraint) 


Lattice QCD 


B K 

fB d \jB Bd 
Am s 


0.88 ±0.08 

(230 ± 12) MeV 
(17.3^;|) ps- 1 


0.90+832 

(231 ± 15) MeV 
(16.3 ±3.4) ps" 1 


0.87 ±0.06 ±0.13 

(230 ± 25 ± 20) MeV 
(15.8 ±2.3 ±3.3) ps" 1 



Table 9: Comparison between the values of the parameters determined in the present analysis 
with those obtained by lattice QCD. In the second row (third row) the results are obtained 
including (not including) the given constraint. 



Values of the parameters given in Table |9| and of the angles reported in Table |^ can be 
considered as reference numbers to which new measurements, using different approaches, 
can be compared in view of identifying new physics contributions. 

Three results will become available in a near future. 

• Analyses of S-meson two body decays have already produced values for 7 which were 
more than 2a larger than the value given in Table ||. Including recent results form 
BaBar and Belle Collaborations, smaller values of the angle 7 are preferred J79|,[50 . 



Before claiming evidence of new physics, a better theoretical understanding of non- 
leptonic decays and more stable measurements are needed. 

Direct measurements of sin2(3 obtained studying J/ipKs events at B factories or 
at TeVatron will provide important tests when the accuracy of these measurements 
will be better than 0.1. 

Finally, a first evidence of B° s — B° s oscillations could still come from LEP and 
SLD Collaborations and accurate measurements are expected from the TeVatron. 
Measured values larger than 25 ps -1 could provide evidence for new physics. 
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